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Abstract 



1 Introduction 



The strong short-range and tensor component in realistic nucleon-nucleon in- 
teractions induces correlations in the nuclear many-body wave functions that 
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cannot be accounted for in an independent-particle approximation as for ex- 
ample adopted in the Hartree-Fock (HF) model. Over the years, various mod- 
els that aim at going beyound the independent-particle model (IPM) wave 
functions, including effects of short-range and tensor correlations, have been 
proposed [1-5]. Despite the strongness of the tensor and short-range force it 
turns out to be a very challenging task to determine the proper measurable 
quantities that can verify the realistic character of those models. Cross sections 
for pion double-charge exchange reactions ^A(7r"'",7r~)^_|_2A have been shown 
to exhibit some sensitivity to dynamical short-range correlations (SRC) but 
cannot yet distinguish between the different model predictions [6,7]. 

Better conditions to study the effects of correlations are predicted for reac- 
tions induced by leptonic probes. In contrast to pion-induced reactions, the 
data are not contaminated by effects related to initial-state interactions and 
the whole nuclear volume, including the interior, is sampled. It is believed that 
with the advent of a new generation of high-duty cycle electron facilities in the 
intermediate energy range (e = 0.5 — 4 GeV) good conditions have been cre- 
ated to explore the dynamics of ground-state correlations with unprecedented 
precision. 

A key function in the study of ground-state correlations in nuclei is the single- 
particle spectral function P{k, E) that gives the probability to remove a nu- 
cleon with momentum k and find the residual A-1 system at an energy E. The 
most dramatic effects of the short-range and tensor correlations on P{k, E) 
are predicted to occur at high momentum and energy [8,9]. One should real- 
ize, however, that these parts of the spectral function belong to the smaller 
probability components in the nuclear wave functions. As a consequence, sig- 
nals of ground-state correlations are likely to produce relatively small cross 
sections. This puts heavy constraints on the experimental requirements when 
performing measurements that aim at probing those correlations. Moreover, 
when exploring the small components in the nuclear wave functions one should 
start worrying about the role of competing mechanisms that could blur their 
effects in the actual cross sections. In this respect it is worth mentioning that 
the frequently quoted relation between the (e,e'p) cross section and the single- 
particle spectral function is based on rather severe assumptions regarding the 
reaction mechanism. Apart from the neglect of any sort of final-state interac- 
tion effect, the initial photoabsorption mechanism is assumed to be dominated 
by one-body operators. Even under quasi-elastic conditions, when this assump- 
tion is believed to have good chances to match reality, there are indications 
for many-body effects contributing to the (e,e'p) cross sections [10-12]. 

Already in the late fifties it was suggested that two-nucleon knockout reactions 
might shed some light on short-range correlations in nuclei. The underlying 
idea is that a photon hitting a strongly correlated nucleon pair will induce two 
nucleons to escape from the target [13]. Two-nucleon knockout reactions have 



been extensively studied using real photons. In these investigations, data have 
been taken in a wide photon energy range (E^ < 1 GeV), covering both the 
proton-proton and proton-neutron knockout channel. Rather than revealing 
information about ground-state correlations, the gross features of the {'y,pp) 
and {'~f,pn) data could be interpreted in terms of meson-exchange and pion 
degrees of freedom [14-17]. Data for electroinduced two-nucleon knockout re- 
actions are rather scarce. Pioneering experiments have been performed at 
NIKHEF-K [18,19]. More data will be collected at MAMI and TJNAF in 
the near future. It is hoped that the additional longitudinal degree of free- 
dom will create optimum circumstances to reveal signatures of ground-state 
correlations. 

One of the principal aims of this paper is to provide a framework in which 
models for ground-state correlations can be confronted with cross sections 
for electroinduced two-nucleon knockout reactions. In our developments we 
leave room for reaction mechanisms that can feed the two-nucleon knockout 
channel and cannot be directly associated with ground-state correlations. The 
major competing reaction mechanism consists of processes in which the pho- 
ton is absorbed on two-, three-, ... nucleon currents. Just as the ground-state 
correlation effects these multi-body currents are a natural manifestation of 
the many-body dynamics of the nuclear system. For the present purposes we 
will restrict ourselves to two-body currents related to one-pion exchange. The 
framework in which those currents are treated, however, is rather general and 
can be extended to include the effects of heavier meson exchange. Apart from 
the exclusive two-nucleon knockout channel, we address also the (e,e'p) reac- 
tion. We concentrate on the high missing-energy region where the effects of 
ground-state correlations are predicted to manifest themselves. The exclusive 
nature of the electroproton reaction at high missing energies cannot be guar- 
anteed. In our theoretical considerations, the starting point will be that the 
single-particle spectral function is not a directly measurable quantity. It will 
be argued that two-nucleon knockout represents an important fraction of the 
continuum (e,e'p) strength and that the (e,e'p) can be considered as a possible 
source of information about many-body dynamics in the nuclear system. 

The outline of this paper is as follows. The specific ingredients of our model 
are introduced in Sec. 2. This includes a discussion of a practical way of cal- 
culating the exclusive (e,e'NN) (Sec. 2.1) and semi-inclusive (e,e'N) (Sec. 2.2) 
cross section. The model assumptions with respect to the meson-exchange and 
isobaric currents are outlined in Sec. 2.3 and those related to ground-state 
correlations in Sec. 2.4. The results of the (e,e'pp), (e,e'pn) and semi-inclusive 
(e,e'p) calculations are discussed in Sec. 3. 



2 Formalism 



2.1 The (e,e'NN) cross section 



In the plane-wave approximation for the incoming and scattered electron waves 
the cross section for triple coincidence reactions of the type (e,e'NaNf,) can be 
written as : 
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(e, e'NaNb) = ——-^kahEaEbfrecCrM 



X 



dEbdVtbdVtade'dVt^, ' ' 4(27r) 

VrWriOa, (pa, Ob, <Ph) + VcWLiOa, (pa, Ob, <Ph) 
+VlWLT{Oa, (pa, Ob, (pb) + VsWrriOa, (pa, Ob, (pb) 

where the Mott cross section (Xm is given by : 



e'^cos^T 



For the cross section (1) we have considered the situation in which the residual 
A — 2 nucleus is created at a fixed excitation energy E^, which is expressed 
relative to its ground-state energy. As a consequence, the integration over the 
energy of one of the escaping particles has been performed. The functions v 
contain all the electron kinematics and read : 



vs = ^. (6) 

The recoil factor in the above cross section reads : 

1 
f = i (7) 



Ea-2 \ ka ka 



where Oab is the angle between the directions of the ejected nucleons. The 
structure functions W are defined in terms of the electromagnetic transition 
operators between the ground state of the target nucleus and the final state 
and are functions of the polar {0) and azimuthal (0) angle of the two escaping 



particles. Our choice for the reference frame and the notation conventions re- 
garding the kinematical variables are summarized in Fig. A.l. In deriving the 
above expression we have not considered any polarization condition for the 
final products and the electrons. Accordingly, the structure functions W in- 
volve a sum over the spin projections {ms^jiris^ and Mf>) of the final products, 
including two escaping nucleons and a residual nucleus created in a specified 
state \ ¥/^-''\e^, JrMr) > : 
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m;P(A = ±1) = {<^f-'>{E^, JrMr)- kafris^; hrris, \ Uq) \ ^o) (12) 
m:P(A = 0) = U>f-^\E^, JrMr); Krris^; kms, \ p{q) \ ^o) ■ (13) 



Here, J\=±i stands for the transverse components of the nuclear current den- 
sity, p is the nuclear charge density and | \l/o > is the ground-state wave func- 
tion of the target system. It goes without saying that computing the above 
matrix elements with a final state characterized by two scattering states can 
be very involving when accounting for the full complexity of the final-state 
interaction (FSI) and that several assumptions have to be made in order to 
keep the calculations feasible. In our model calculations we have put special 
emphasis on the orthogonality condition between the initial and final states 
thus avoiding spurious contributions entering the matrix elements. This is of 
particular importance in view of the fact that on several occasions we will 
make integrations over relatively large parts of the available phase space for 
the final products. In two-nucleon knockout from finite nuclei the final-state 
interaction involves apart from mutual interactions between the escaping nu- 



cleons also distorting effects from the residual A — 2 system. Another aspect 
of the FSI that deserves further attention is the role of multi-step processes. 
An obvious way to take these into account is a coupled-channel calculation. 
In the two-nucleon emission case, however, lots of channels are expected to 
contribute which makes a coupled-channel calculation very challenging. In the 
absence of a consistent and handleable theory to describe the full complexity 
of the FSI, we adopt the view that rather than accounting for only part of 
the FSI it is maybe better to work in a direct knockout model. Or put in 
other words, we assume that the two nucleons involved in the photoabsorp- 
tion mechanism will escape from the target nucleus without being subject to 
inelastic collisions with the core. This is the basic assumption of the so-called 
"spectator approximation" . In the shell-model picture, direct emission of two 
nucleons will leave the residual nucleus in a two-hole (2h) state relative to 
the ground state of the target nucleus. Even in a direct reaction model, the 
residual A — 2 nucleons will distort the wave function of both escaping nucle- 
ons. In our model calculations this FSI effect is implemented along the lines 
explained in Ref . [20] . There it is shown that a proper A-body wave function 
with two asymptotically free nucleons and a residual nucleus generated in a 2h 
state can be reached by performing a partial wave expansion in terms of the 
two-hole-two-particle eigenstates of a mean-field one-body Hamiltonian. For 
the sake of completeness we mention the obtained expression for the required 
type of final wave function : 
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xF4^(l],)F;„^,(a) < Imi-nis^ \ jm >< jmj'm' \ JiMi > 

X < I'mii-rus^ I j'm >< JrMrJiMi \ JM > 

X I {hh'Y'EJn- {p{eJj)p'{e,l'f))J^; JM > . (14) 

The model that we adopt to account for the FSI is thus based on a partial wave 
expansion in terms of two-hole-two-particle (2h2p) eigenstates of a mean-field 
potential. These wavefunctions are defined according to 
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The continuum or particle {p) eigenstates of the mean-field potential are char- 
acterized by p{elj). The energy at which the partial waves p{elj) are calculated 



is determined by the momentum of the emitted nucleon : e^ = k'^ /{2Mn). The 
central and Coulomb phase shifts are denoted by 5i and o";. Throughout this 
work we consider isospin not to be a good quantum number, in the sense 
that at the level of the single-particle states we discriminate between protons 
and neutrons. The single-particle wave functions are constructed through a 
Hartree-Fock calculation with an effective Skyrme type interaction [21]. For 
the protons a Coulomb part is added to the mean-field potential. 

The state in which the final state is created is determined by | {hh')~^ExJR >, 
with Ex the excitation energy with respect to the ground-state energy of the 
residual nucleus. A schematic drawing that illustrates the basic idea behind the 
expansion (14) is shown in Fig. A. 2. It should be noted that in reality only a 
fraction of the final-state wave function will be of two-hole nature. Indeed, the 
many-body character of the residual system will make the two-hole strength 
to be fragmented over a wide energy range. In order to account for this nuclear 
structure effect, "spectroscopic factors" have to be considered when calculating 
cross sections for decay to specific states. The nuclear-structure aspects which 
are expected to play a role in the ^^0(e,e'pp) reaction have been the subject 
of a recent study reported in Ref. 



At a first sight, the usefulness of Eq. (14) is questionable as all expansions 
extend to infinity. When performing the calculations in coordinate space, finite 
expressions for the transition matrix elements can be obtained by performing 
an angular decomposition of the transition operators {J\{q), p{q))- This leads 
to the well-known expansion of the nuclear four-current operator in terms of 
the Coulomb (Mj"^'), electric (Tfj^) and magnetic {Tj'^^) multipole operators 
(J= V27n) 



J±,(g) = -v^ E r^J [Tf^M ± T-f{q)) (16) 

j>i 

p{q) = V^Y.''jMjT\^)- (17) 

J>0 

After inserting the expansions for the final state (14) and the nuclear charge- 
current operator (16), the transition matrix elements of Eqs. (12) and (13) 
can be cast in the closed form : 



<(A = ±i) = -v^E^''^ E E E 

J>1 Imijml'miij'm' JiMi 

x(47r)2(-2)'+'' ^^e-^(^'+'^'+''"+'^") 

1 1 



X < jmj'm' I JiMi > ^ —. < Jr- MrJX \ J^Mi > 

X [-M J';M'(-/i, J. Jr) + ^M;ZA.h, J, Jr)] 



m 



{x=o) = v^j:^'jt. e e 



J>0 ImijmVmiij'm' J\Mi 

^ -j(<5;+a;+5,,+a(,) 



2MN^/kah 
xYirai{yta)Yi>mi,{S^b) < Imi-rUs^ \ jlTl X I'mv-iTis, \ fm' > 

X < jmj'm' I JiMi > ^ ^—. < Jr- MrJO \ JiMi > 

Ji 
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where the matrix elements Ai have been defined according to : 

-MS'7/f'(-/i, J. Jr) =< p{^i.ij)p\^J'f); J^Wrf'-^^^iq^hh'; Jn > 

_^_iy.+J.'+Jn < p(e,lj)p'{ej'j'); J4Tf'"'''<^{q)\\h'h; Jr > , (20) 

M;";;,^!, J, Jr) =< PiebWieJ'f); J4MrKq)\\hh'; Jr > 
_^_iy.+Jn'+jR < p{e,lj)p'iej'j'); J,\\Mr\q)\\h'h; Jr > . (21) 

It should be stressed that the above expressions are general and that no as- 
sumption whatsoever regarding the nature of the inital photoabsorption mech- 
anisms has been adopted. The contributing photoabsorption mechanisms will 
be discussed in Sects. 2.3 and 2.4. 



2.2 The semi- exclusive (e,e'N) cross section 



In comparison with exclusive (e,e'NN) reactions the semi-exclusive (e,e'p) 
channel is more attractive from the experimental point of view. For model cal- 
culations, however, the semi-exclusive channel can be more challenging than 
exclusive 2N knockout calculations. Indeed two-nucleon knockout is expected 
to feed the semi-exclusive (e,e'p) channel as soon as the thresholds are crossed. 
The calculation of the two-nucleon knockout contribution to the semi-exclusive 
channel involves integration over the phase space of the undetected escaping 
nucleon (either a proton or a neutron). In what follows we will exploit the fact 
that most of the derivations for the two-nucleon knockout cross sections in 
previous sections are done in coordinate space to perform some of these inte- 
grations analytically. In this way we facilitate the calculation of semi-exclusive 



cross sections enormously and develop a framework that can be further ex- 
ploited to investigate the role of the dynamical nucleon-nucleon correlations in 
the semi-exclusive channel. As most of the dynamical correlations are believed 
to be of two-body nature, one could expect that they manifest themselves most 
clearly in the two-nucleon emission channel. It is a major challenge to trace 
those kinematical regions where the semi-exclusive channel is predominantly 
fed by emission of two nucleons. Such investigations start from considering the 
contribution from (e,e'pp) and (e,e'pn) to the semi-exclusive channel : 
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An obvious way of determining the right-hand side of the above expression is 
calculating the 2N knockout cross sections for a grid of kinematical conditions 
that covers the full phase spaces dVtp'dEpi ((e,e'pp) contribution) and dVtndEn 
((e,e'pn) contribution). Those cross sections can then be integrated to obtain 
the semi-exclusive strength that is attributed to 2N knockout. Such a proce- 
dure would be cumbersome and would consume a lot a computing time. We 
present an alternative method that avoids those numerical integrations at a 
very small cost for the accuracy of the calculated cross sections. We start with 



remarking that most of the two-nucleon emission strength tet-jo jip jo j /jo 
(a=p',n) is confined to a relatively small fraction of the complete dfladfip phase 
space : ever since the original work by K. Gottfried [13] it has been clear that 
2N knockout reactions predominantly occur in back-to-back situations (quasi- 
deuteron kinematics). Consequently, for a particular dEpdVtp (semi-exclusive 
situation) most of the (e,e' N^ p) strength will reside in a restricted angular 
range for the second nucleon N^ that remains undetected. In this restricted 
dVta area the value of the undetected nucleon momentum, which is determined 
from energy-momentum conservation, will vary very slowly so that to a good 
approximation it can be replaced by an "average" value k'^"^ . This average 
momentum can be determined by imposing quasi-deuteron kinematics (which 
is equivalent with considering the situation with zero recoil momentum ) : 

q-kr-kp = ^, (23) 



where kp is the momentum of the detected proton. After introducing this 
"average momentum", the integration over dVtpi and dVtn in Eq. (22) can be 
performed analytically. Indeed, after inserting the expressions (18) and (19) 
into the structure functions W{9a,(f>a,(^b,(pb) of Eqs. (8-11) it can be shown 
that : 
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where the function B was defined according to : 
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These expressions allow to compute directly the semi-exclusive cross sections 
as a function of the kinematical variables {6p,(f)p,Tp) without explicitly calcu- 
lating the intermediate (e,e'NN) angular cross sections. In the actual calcu- 
lations we account for the spreading of the two-hole strength in the energy 
spectrum of the A-2 system. This is done with the aid of a model explained 
in Ref. [23]. 



2.3 Meson-exchange and isobaric currents 



In an independent-particle picture, the most direct source of two-nucleon 
knockout strength would be the coupling of the electromagnetic field to the 
hadronic two-body currents. These currents are mediated by the mesons which 
lie at the basis of the one-boson exchange picture of the nucleon-nucleon inter- 
action. Here, we do not make any attempt to consider medium-corrections for 
the meson fields. In constructing the two-body currents, our starting point is 
the one-boson exchange (OBE) picture of the nucleon-nucleon interaction as 
it has been derived from nucleon-nucleon scattering data. The long-range part 
of the OBE potential is well established and dominated by pion exchange. In 
coordinate space it takes on the well-known form as a sum of a spin-spin and 
tensor interaction term [24] 
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6 \ rriT, 
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The current operator satisfying the continuity equation with this potential 
can be constructed from general field-theoretical methods. We have consid- 
ered pseudovector nNN couphng and treated the two-body currents in the 
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non-relativistic limit [25]. This results in the well-known 7r-exchange current 
operators : 
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where the first two terms refer to the so-called Seagull current and the last 
term, which is quadratic in the pion propagators, to the pion-in-flight current. 



In addition to the meson-exchange currents (MEC) derived from the OBE 
potential we have also considered processes in which nucleon excitations A33 
are created (Fig. A. 3). Here, the model dependency is intrinsically larger as 
the (transverse) currents cannot be constrained by the continuity equation. In 
order to construct the currents associated with A33 creation with subsequent 
pion exchange, we also rely on field-theoretical methods and adopt a non- 
relativistic approach. The ttNA coupling is considered in the standard form : 
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where S and T denote the spin and isospin | — ^ | transition operators. In the 
7A^A coupling we retain solely the magnetic term 
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with A the external electromagnetic field and F^{ql) the electromagnetic form 
factor of the delta for which we used the standard dipole form. The (marginal) 
electric quadrupole term in the 7A^A coupling is neglected throughout this 
work. With the above coupling lagrangians we arrive at the following expres- 
sion for the A33-current with 7r-exchange : 



f r;7 -7 -7 ^ - ^ f^NAfnNNfnNA jp / 2^ 

JwA[q,qi,q2) - ^ ^A[qnj 
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non—res 
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02^ + ml 
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Qi—^ 
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X g 



(32) 



The vrA^A and 7A^A coupling constants are taken from ref. [26] : J'^t^na I 
4 TT = 0.37 , /^AfA = 0.12. In the above expression for the isobaric current 
G^* denotes the A33 propagator for the resonant diagrams (Fig.A.3(a) and 
(c)). The corresponding propagator for the so-called non-resonant diagrams 
(Fig. A. 3(b) and (d)) is denoted by Q'^^-'^^^ _ As we will consider energy trans- 
fers in the resonance region, special attention has been paid to constructing 
the A33 propagators. A free A33 excitation obtains a width through ttA^ decay. 
The corresponding propagator would then read 



Gres 
A 



-E2' + Ma- f r^"' 



(33) 



where E^^^ is the intrinsically available energy for the resonance and the width 



r^^^ becomes [26,27] 



T^res 
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(34) 



where j^ is the decay momentum in the center-of-mass (c.o.m) frame of the 
TiN system and ^/s the total c.o.m. energy of the pion and nucleon. In the 
medium, however, the nN decay will be blocked by the Pauli principle. Various 
pion-nucleus [28], photoabsorption [29] and inclusive (e,e') [30,31] experiments 
have pointed towards other strong medium modifications of the A33 resonance. 
Indeed, the inclusive A(e,e') spectra show a pronounced broadening and damp- 
ing of the resonance in comparison with A times the free electron nucleon cross 
section. It is common belief that this is mainly due to the coupling of the dom- 
inant A -^ ttN decay with the AA^ -^ NN channel. It is precisely the latter 
channel which is under investigation in two-nucleon knockout experiments. 
The A-hole model [30] has been applied successfully to the description of in- 
clusive (e,e') reactions in the A33 resonance region [31]. In the A-hole model 
one adopts a dynamical description of isobar propagation in the nucleus. This 
results in a propagator of the form [30] 



A 



rpres 



+ M^ 



t 'pres 
2^ A 



5W + W^ + W, 



(35) 
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where SW accounts for Pauli blocking and W-^ for coupling of the resonance 
to tt" and the nuclear ground state. These two terms tend to cancel each other 
so that Wsp will induce the major correction with respect to the free delta 
propagator of Eq. (33). The Wgp is a semi-phenomenological parametrization 
of the coupling of the delta with the more complicated channels. In the orginal 
formulation of the delta-hole model the Wgp was chosen to have a density- 
dependent central and spin-orbit part. Chen and Lee [32] have shown that an 
equally good description of the ^^C(e,e') cross sections could be obtained by 
simply assuming that 

Wsp[MeV] = -30 - 40 i . (36) 



This prescription is in agreement with the results from recent total photoab- 
sorption measurements [29] from which it was concluded that the A mass 
and width increases with growing nuclear density. The values for the A mass 
and width for ^^C quoted in Ref. [29] are in agreement with the above pre- 
scription. In what follows we will adopt the procedure (36) to account for the 
medium effects on the A33 propagation. Recently, this approach has been ap- 
plied with some success to the description of the photonenergy dependence of 
the ^^C(7,pp) and (7,pn) cross sections in the resonance region [33]. 

It is well known that the delta peak in the inclusive (e,e') cross section appears 
at a higher energy transfer u than 300 MeV which is typical for real photon 
induced reactions. It is convenient to introduce the equivalent photon energy 
K to produce the same nuclear excitation energy as the virtual photon (g, u) 



A real photon with energy K produces the same ttN c.o.m. energy as a virtual 
photon (g, u) . 

As we are dealing with off-shell nucleons the photon energy is not completely 
available for internal excitation of the nucleon. As pointed out in Refs. [34,35] 
a reasonable substitution for E^"^ is : 

(EX")' = (Mn - ehf + 2K{Mn - en) , (38) 



where th is the average binding energy of the nucleon on which the pion is 
reabsorbed. Remark that e^ depends on the shell in which the nucleon is 
residing. 
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For the non-resonant diagrams (Fig. A. 3(c) and (d)) the propagator becomes 

rinon—res ^'^Q'l 

^ - _^non-res ^ ^^ , l-^yj 



with 



(T7inon—res\ 
Ea ) 



\ 




(40) 



In the static hmit (small uj and q), one obtains : 

1 

Ma- Mn ' 



/^res /~mon—res I AT\ 



in which case the isobaric current operator of Eq. (32) reduces to the static 
operator form as e.g. derived in Ref. [36] . By no means the static limit should be 
considered realistic as soon as one is approaching equivalent photon energies 
K that probe the real resonance region. 

In fitting nucleon-nucleon scattering data in terms of a particular One-Boson 
exchange potential it is a common procedure to regularize the tiNN vertices 
for the finite size of the hadrons and the complexity of physical mechanisms 
that are thought to happen at distances which are short as compared with the 
average range of the pion. This is frequently done by introducing a monopole 
hadronic form factor. 

-^ (42) 



at each ttNN vertex. For the actual calculations, these hadronic form factors 
have also been introduced in the meson-exchange and isobaric current oper- 
ators of Eqs. (29) and (32). Throughout this work, a pion cut-off mass A^ of 
1250 MeV has been used [24]. 



2.4 Ground-state correlations 



In all of the above considerations we have adopted the independent particle 
model (IPM) in which all nuclear wave functions are cast in a Slater deter- 
minant form. In this Subsection we present a method to impose corrections 
to that picture. To that purpose we rely on a perturbation expansion method 
[37,38] to calculate transition matrix elements between many-particle wave 
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functions that have been corrected for correlation effects that go beyond the 
IPM. Here, correlations in the nuclear wave functions are implemented via 
a technique inspired by the so-called correlated basis function (CBF) theory 
[39] in which correlated wave functions \E' are derived from their (uncorrelated) 
independent particle limit through the operation of an operator Q. The latter 
corrects the Slater determinant \E' for short-range and other correlations not 
accounted for in the IPM : 

— ^1 ^) 

^) = I = ■ (43) 

'^ I g^g I ^' 



In the CBF theory, the operator g takes on the form of a symmetrized product 
of two-body correlation operators : 



g = s 



A 



.i<j=l p=l,6 



(44) 



where fij = fi — fj, Rij = {fi + fj)/2 and S is the symmetrization operator. 
The Of- are the isoscalar and isovector isospin operators of the scalar, spin 
and tensor type [40] : 



df^e{l{p = l),r, ■ T,{p = 2), a, ■ a,{p = 3), (45) 

(Tj ■ ajfi ■ fj{p = 4), Sij{p = 5), SijTi ■ fj{p = 6)} . 

The different components in the operator g reflect the fact that the nucleon- 
nucleon force is a function of the spin and isospin orientation of the interacting 
particles. Of all of the above components the central {p = 1) and tensor {p = 6) 
operator have been notified [5,41] to induce the largest correlation corrections 
to the nuclear Slater determinants obtained in the IPM approach. 

Throughout this paper we will restrict ourselves to the effects induced by the 
central operator {p = 1), that induces short-range correlations (SRC) to the 
nuclear wave functions obtained in an IPM approach. It should be stressed, 
however, that the techniques expanded in this paper can also be applied to the 
other terms. In most of the calculations that address the correlation function 
fij~ {"^ijjRij) it is assumed that the functional dependence can be restricted 
to the relative coordinate of the interacting nucleon pair Vij =\ fi — fj \. This 
assumption, which considerably simplifies the calculation of the correlation 
functions, can be justified by considering that the range of the correlation 
function jf~ is small with respect to the surface thickness of the nucleus. 
Accordingly, the dependence of ff~ on the c.o.m. coordinate Rij is anticipated 
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to be small. As a result of retaining solely the corrections induced by the short- 
range correlations, the correlated nuclear wave functions used throughout this 
paper read 

^ix^,....,XA)= n /r'M^(^i---^A)/v^' (46) 

where x is a shorthand notation for the radial and spin coordinates and A^ is 
the normalization factor 



dxi....dxA'^\xi,....,XA)i n fl 



Within the adopted assumptions an arbitrary transition matrix element be- 
tween a correlated ground state | ^o > and a final state | \1'/ > can be 
rewritten as 

^; I O I ^o) = -y^ (^f I O'" I ^o) , (48) 



with 



0^^^= n (l-^(r.,))^0 n (l-^M) , (49) 

i<j=l k<l=\ 



where we have introduced the central correlation function g which is defined 
according to g{rij) = 1 — ff~ {vij). In the absence of short-range corrections 
the correlation function g{rij) would simply be zero. The effective operator 
approach as formally written in Eq. (48) has the marked advantage that we 
can rely on standard techniques, like second quantization, when calculating the 
transition matrix elements between correlated many-body states. The adopted 
approach also allows to put the short-range correlations on the same footing 
as other mechanisms that can feed the electroinduced two-nucleon knockout 
channel. In this way it will become easier to evaluate the sensitivity of the 
calcaluted cross sections to the SRC relative to the other mechanisms. The ef- 
fective operator of Eq. (49) is an A-body operator which makes the exact calcu- 
lation of "correlated" transition matrix elements of the type (48) not feasible. 
Various cluster expansions, however, have been developed to approximate the 
matrix elements of Eq. (48) [38,42]. These techniques usually address matrix 
elements of the type < \I'o | O | \E'o >, where ^o is the correlated ground state 
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of the nuclear system. Here, we are facing a quite different situation in the 
sense that we are deahng with transition matrix elements between wave func- 
tions of totally different origin, namely \l/o and a final state with two nucleons 
residing in the continuum. Given that two-nucleon knockout calculations are 
already quite cumbersome in a pure IPM, there is a practical limitation on 
the number of terms in the cluster expansion that can be accounted for. This 
does not necessarily mean, however, that the main effects of the SRC cannot 
be calculated to a high degree of accuracy. Most of the cluster expansion tech- 
niques developed for the calculation of the "correlated" matrix elements rely 
on an expansion into the different orders of the central correlation function g. 
Here, we will restrict ourselves to the terms linear in the correlation function 
g. Higher order terms require three and more correlated nucleons. Given that 
the central correlation effect is very short ranged, this type of multi-nucleon 
correlations could be expected to be less important than two-body correlations 
at normal nuclear densities. In all forthcoming derivations we assume that the 
transition operator O in expression (48) has a one- and two-body part : 



d = Y.o^%) + Y.o^%, 



3) 



(50) 



i<j 



Within the adopted assumptions the effective transition operator can then be 
written as : 



ij) 



Oe// = ^OW(^)+EOP](^,j)-E 0W«+0W(J) 9{r. 

i i<j i<j '- 

- E [oW(0^M+OW(j>M+OW(A:)^(r.,) 

i<j<k 
i<j 

i<j<k 

+d^'\t, k)g{r,k) + 0[2](j, k)g{n,) + 6"^%, k)g{n,) 

- E W\h3)9{rH)+d^%,k)g{r,i) + m{i,l)g{r,,) 

i<j<k<l'- 

+d^^\j,k)giru) + d^'\j,l)gir,,) + d^'\k,l)gin,) 



(51) 



For the sake of brevity we did not write the terms in g"^, that refer to final- 
state correlations, in the above expression. Formally these contributions have 
exactly the same form as the terms generated by the initial-state correlations. 
The above expression clearly illustrates how the short-range effects are put 
on the same footing as the other contributions in the photoabsorption mech- 
anism. The effective transition operator has apart from the terms occuring 



in the "original" operator O additional terms that depend on the correlation 
function g. The implications of those terms for the two-nucleon knockout re- 
action have been illustrated in Fig. A. 4. In the absence of correlations that 
go beyond the IPM only diagrams of the type (a) would contribute to a di- 
rect two-nucleon emission process. Diagram (a) corresponds with a genuine 
two-body absorption process as e.g. encountered in the situation in which a 
photon is coupled to a charged meson exchanged between two nucleons mov- 
ing in a mean-field orbital. Diagram (d) is of the same form as diagram (a) 
and can be interpreted as a short-range correction to the two-body current 
contributions from diagram (a). It will introduce a decreased probability at 
short internucleon distances for two-body photoabsorption processes to occur. 
As those short-range events are already heavily cut by the introduction of 
hadronic form factors they have a rather small impact on the calculated cross 
sections [43]. 

Diagrams (c),(e) and (f) are typical events that are linear in the correlation 
function g and involve three and four-body operators. Given the numerical 
complexity of two-nucleon knockout calculations we will restrict ourselves to 
the two-body operators that arise from diagrams of the type (a) and (b). This 
is often referred to as the Single-Pair Approximation (SPA). Within the SPA, 
one retains those terms from the above expansion that refer to processes in 
which the photon hits two correlated nucleons as a result of which they are 
emitted from the target nucleus. In the calculations we do not account for 
reactions in which the photon couples to a nucleon as a result of which one 
(or two) nucleons in the recoiling (A-1) system are emitted. Moreover, we 
did not consider the final state correlations. Indeed, within the context of the 
SPA inclusion of the final state correlations would mean that one considers 
the correlation between two nucleons residing in a scattering state. This type 
of correlation could be expected to be small. To the best of our knowledge, 
no calculations are available that address the correlation function for the con- 
tinuum states appearing in the final state of the reactions under study. Most 
investigations into the correlation function have solely addressed the nuclear 
ground state. 

It is worth stressing that in the absence of ground-state correlations the one- 
body photoabsorption operator O'^^ would not contribute to the direct din- 
ucleon emission process. It is only after introducing the central (or Jastrow) 
correlations that it starts contributing and can be formally treated like e.g. the 
MEC Referring to Eq. (51) it seems technically possible to write down effec- 
tive many-body currents that arise from combining one-body photoabsorption 
and ground-state correlation effects. In the SPA, the effective two-body current 
that accounts for the short-range corrections in the ground state reads 

jfRcihj) = - [J^\i) + J^%)) 9{n,) . (52) 
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For the one-body current jt^l we consider the standard form of the Impulse 
Approximation (lA) that has a convection and magnetization part. Conse- 
quently, 



^SRCV'^h 



2iM, 



N 



VJ 



^3 



2iMN 
'2Mn 



2M, 



N 



WiSif- fi) + 5(f - fj)ViJ g{r, 
Vj5{f- fj) + 5{f- rj)Vj\ g{rij) 

5{f-fi)V X aig{rij) 

S{f-fj)V X ajg{rij) . 



Similarly, the effective two-body charge density becomes : 

PsRC^hJ) = - {ei6{f- fi) + ej6{f- fj))g{rij) 



(53) 



(54) 



In order to facilitate the calculation of the matrix elements, the central corre- 
lation function is expanded according to : 



g{ri2) = Y. 



Im 



At: 
2/ + 1 



gi{n,T2)Y,Un,)Yim{^2) . 



(55) 



The different partial-wave components can be directly obtained from : 

+1 



9i{n,r2) 



2/ + 1 



dxPi(x)g i yrf + r^ — 2rir2X J . 



(56) 



The expressions for the two-body transition matrix elements (20) and (21) 
with the effective two-body operators of Eqs. (53) and (54) are given in Ap- 
pendix A. These matrix elements represent are at the basis of our calculation 
of the SRC contribution to exclusive (e,e'NN) and semi-exclusive (e,e'N) cross 
sections. Results of those calculations will be presented in the forthcoming sec- 
tions. In the calculations which will be presented below we have used several 
forms of the central correlation function g. Theories that start from princi- 
pal grounds, seem to produce correlation functions that are rather soft. In 
this context we refer to the Fermi Hypernetted Chain (FHNC) calculations 
of Refs. [3] and [44], the Monte Carlo calculations that were e.g. reported in 
Refs. [2] and [45] and the recent calculations performed in a translationally in- 
variant framework [5]. The object of these studies is usually the ground-state 
properties of finite nuclear systems. Most of these calculations include also 
higher order cluster terms. 
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To set some kind of upper limit on the effect of tlie SRC for the calculated cross 
sections, we will also present results with a pronounced hard-core correlation 
function that is due to Ohmura, Morita and Yamada (OMY) [46] : 

.0 r < 0.6 fm 

fir) = { - ,(57) 

[1 — exp {—^i^{T — cY)] [1 + 7exp (yU^(r — c)^)] r > O.Qfm 



with fi = 1.118 fm"^ and 7=2.078. We do not consider this as a realistic 
correlation function but use this function to maximize the effect of the short- 
range correlations relative to the meson-exchange and A degrees of freedom. 



3 Results and discussion 



3.1 Exclusive (e,e'pp) and (e,e'pn) reactions 



The most elementary (e,e'NN) cross sections (c.s.) that can be calculated are 
those for excitation of the residual nucleus in a specific (discrete) state. Even 
with fixed the electron kinematics, one is left with an angular cross section 
depending on five independent variables. An obvious choice for the latter are 
the two solid angles and one of the kinetic energies of the escaping particles : 
{fla,^b,Ta). In what follows we will restrict ourselves to in-plane kinematics 
which fixes the azimuthal angle of the two detected nucleons at 0° and 180° 
degrees. In Figs. A. 5 and A. 6 we show -'^^C(e,e'pp) and -'^^C(e,e'pn) angular 
cross sections for creation of the final nucleus in a state with | (1^3/2)"^, Jr > 
two- hole structure and Jr = O"*" and 2"*". The considered kinematical conditions 
are typical ones for the 800 MeV AMPS (Amsterdam) and MAMI (Mainz) 
electron accelerators. In order to further reduce the number of independent 
variables, for the calculations presented here we have considered the situation 
in which both escaping nucleons have equal kinetic energies at zero recoil 
momentum Pa-2 

Pa-2 = q- ka- h . (58) 



This means that the kinetic energy for one of the escaping particles was fixed 
at T = (cj — S'2Ar)/2, with S2N the threshold energy for 2N emission. With 
this choice we are left with two independent variables, the polar angles of the 
escaping nucleons, against which the angular cross sections can be studied. In 
Figs. A. 5 and A. 6 the polar angles vary between 0" and 360°. Within this 
convention a typical back-to-back situation in the LAB frame would be e.g. 
{6a = 40°, 6'b = 220°). The thresholds for proton-proton and proton- neutron 
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knockout from ^^C are respectively 27.2 and 27.4 MeV. Obviously, the two- 
nucleon knockout strength is residing in a very small part of the phase space, 
namely those situations in which both nucleons escape in an almost back-to- 
back situation. This type of behaviour can be easily explained by referring to 
the factorized model for two-nucleon knockout [13,47]. In such an approach it 
can be shown that the angular cross sections are proportional to the c.o.m. 
distribution F{P) of the active nucleon pair. As the F{P) is a sharply decreas- 
ing function with c.o.m. momentum \ P \ {P = —'pA-2), two-nucleon knockout 
will preferentially occur in those situations in which the momentum P is kept 
small. From energy-momentum conservation arguments, this corresponds with 
back-to-back emission. Despite the fact that in the unfactorized model some of 
the approximations which are at the basis of the factorization are not made, we 
still find the c.o.m. distribution of the pair to be a driving mechanism behind 
two-nucleon knockout. Besides this back-to-back dominance, it is clear that 
other variables do play a role. Inspecting Figs. A. 5 and A. 6 it is obvious that 
the angular momentum of the residual nucleus has a large impact on the shape 
of the angular cross section. The angular distributions for pp and pn knockout 
leading to a 2"'" state are much wider than those for the corresponding 0"*" state. 
Moreover, also the details of the leading photoabsorption mechanisms reflect 
themselves in the shapes of the angular cross sections. Even for creation of the 
final nucleus in a particular J^ and corresponding kinematical conditions, the 
proton-proton and proton-neutron cross section show some differences which 
can be related to the reaction mechanisms contributing to the cross sections. 
In this context, we want to remind the reader of the fact that the meson- 
exchange (Eq. (29)) and isobaric (Eq. (32)) currents have a peculiar isospin 
dependency which makes them behaving differently in the proton-proton and 
proton-neutron channel. 

One of the main tasks of this paper is to investigate the sensitivity of the 
(e,e'NN) cross sections to dynamical short-range correlations. For that purpose 
we have calculated ^^0(e,e'pp) and (e,e'pn) angular cross section for emission 
out of the p-shell orbitals. A measure of the relative importance of the ground- 
state correlations in the whole reaction process, is the contribution from the 
longitudinal channel to the differential cross sections. Indeed, into first order 
the longitudinal channel is free from isobaric and meson-exchange current 
contributions. In the presented model calculations the longitudinal strength is 
exclusively the result of the existence of dynamical short-range correlations. 
A measure for the degree of longitudinal polarization of the virtual photon is 
the parameter e^ 



Vt 




(59) 
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In Figs. A. 7- A. 10 we have plotted some corresponding proton-proton and 
proton-neutron knockout angular cross sections at two values of the energy 
transfer. These have been chosen as to probe a kinematical regime well below 
and above the A resonance region (the equivalent photon energies K are 198 
and 352 MeV respectively). The large incident electron energy (e=1.2 GeV) 
and small electron scattering angle {6^=12 °) considered make the longitudinal 
polarization ei of the virtual photon large. For the results of Figs. A.7-A.10, 
we have considered the situation in which the residual nuclei ^^C and ^^N are 
created in a | (lp3/2)~^(lpi/2)~"'^; <//? = 2+ > state. According to the nuclear- 
structure calculations of Ref. [22] this type of configuration would be prefer- 
entially populated in ^^0(e,e'pp) reactions. Furthermore, we have considered 
the situation in which both escaping nucleons have the same kinetic energy at 
I 'Pa-2 1= 0. The proton- neutron channel is obviously a much stronger channel 
than the proton-proton emission channel. Inspecting Figs. A. 8 and A. 10 it 
is clear that proton-neutron emission is dominated by the transverse channel 
and that central ground-state correlations play only a marginal role. 

The weaker proton-proton emission channel exhibits a stronger sensitivity to 
the ground-state correlations. Over the whole, however, the isobaric currents 
dominate the angular cross sections. There are some parts of the phase space, 
however, for which the longitudinal cross section is a sizeable contribution in 
the full cross section, reflecting the sensitivity of the reaction to ground-state 
correlation effects. The latter are e.g. playing a major role for kinematical 
conditions that correspond with one of the escaping nucleons moving along 
the direction of the three-momentum transfer q. 

All of the previous angular (e,e'NN) cross sections have been obtained with 
the central FHNC correlation function. The latter is of the Gaussian type 
5^(^12) = a;e~^''i2. For ^^O, the parameters were calculated to be q;=0.53 and 
/3=1.52 fm~^ [3]. For our purposes, the parametrization for ^^C was assumed 
to be identical to the one for ^^O. 

We now aim at studying the sensitivity of the calculated (e,e'pp) cross sections 
to the choice of the Jastrow correlation function. To this end, we do not longer 
consider the situation in which the residual nucleus is created in a particular 
Jji but plot the cross section which is obtained after (incoherently) adding 
all possible final angular momenta Jr. The number of situations which have 
to be considered depends on the quantum numbers {IhjhJh'jh') of the orbits 
from which the nucleons are escaping. For example, for proton-proton emis- 
sion out of the (lp3/2, lpi/2) orbits the plotted cross section is obtained after 
adding the contribution from Jr = 1+ and 2"'". For the results of Fig. A. 12 we 
consider the so-called coplanar and symmetrical situation. This type of kine- 
matics has been schematically sketched in Fig. A. 11 and corresponds with the 
situation in which both nucleons are escaping with equal kinetic energy and 
opening angle with respect to the direction of the transferred momentum. In 
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Fig. A. 11 we show also the typical behaviour of the c.o.m. momentum P in 
such kinematics. So given the dominance of the c.o.m. distribution in pro- 
ducing the general features of the two-nucleon knockout cross sections, one 
would expect the angular cross section to peak around an opening angle of 
70 degrees and fall off quickly to smaller and larger opening angles. This is 
precisely what is observed in Fig. A. 12. Another observation is that with the 
central correlation function of Ref. [2], which is obtained with Monte-Carlo 
techniques, the Jastrow correlations are hardly affecting the differential cross 
sections that are obtained when solely accounting for the A degrees of free- 
dom. In the remainder of the paper, the correlation function of Ref. [2] will 
be referred to as "MC" . A more pronounced effect is found with the Gaussian 
FHNC correlation function. The short-range contribution is further noticed 
to depend on the single-particle levels from which the protons are escaping. 
Whereas the Gaussian FHNC correlation function slightly increases the cross 
sections for (lpl/2)^^ and (lj93/2)^^ knockout, the effect is quite sizeable for 
2N knockout out of the mixed {lpl/2)'^{lp3/2)'^ orbits. The most spectac- 
ular shell dependence is noticed for the hard-core OMY correlation function. 
Counting the number of possible proton-proton pairs one would expect the 
absolute cross sections for the different two-hole configurations of Fig. A. 12 to 
be like 6 {{lp3/2y^y.8 {{lp3/2)~\lpl/2y^y.l{{lpl/2y^). It is obvious that 
the number of pairs can only serve as a rough guide to predict the relative 
amount of strength that goes into the different two-hole combinations. 



In Fig. A. 13 the four terms (longitudinal, transverse, longitudinal-transverse 
and transverse-transverse) contributing to the ^^0(e,e'pp)((lp3/2)~-'^(lpl/2)^-'^) 
angular cross sections of Fig. A. 12 are shown for different choices of the central 
correlation function. It is clear that the absolute magnitude of the cross section 
is dominated by the transverse (T) and transverse-tranverse (TT) terms that 
show very little sensitivity to the central correlation effects. The longitudinal 
(L) and longitudinal-transverse (LT) structure functions, on the other hand, 
exhibit a very strong sensitivity to the choice of the correlation function. The 
hard-core OMY correlation function produces longitudinal (e,e'pp) strength 
that is typically one order of magnitude larger than the strength generated by 
the (realistic) soft-core MC and FHNC correlation functions. Remark further 
that the MC and FHNC correlation functions produce remarkably different L 
and LT structure functions. This is a rather surprising result as their func- 
tional dependence on the relative coordinate of the two interacting particles 
is rather similar. 
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3.2 Semi- exclusive (e,e'p) 



In the plane-wave impulse approximation (PWIA) the (e,e'p) cross section can 
be cast in the form 

EpaepP{\Pra\,Era) , (60) 



de'dVte'dVtpdEp 



where Ep (pp) is the energy (momentum) of the detected proton and a^p the el- 
ementary cross section for electron scattering off an off-shell proton. It should 
be noted that a similar type of factorization can be pursued for the (7,p) 
reaction. The spectral function P(|pm| , -E'm) is related to the probability of 
removing a nucleon with momentum pm from the target nucleus and finding 
the residual nucleus at a missing energy Em- The spectral function is of funda- 
mental importance for the understanding of the many-body dynamics of finite 
nuclei and has been the subject of many theoretical considerations [42,48,49]. 
Eventhough a wide variety of spectral functions is available in literature, most 
calculations point towards the IPM providing a fair description of P(A;, E) at 
low k and E. As far as the spectral function is concerned, the deviations from 
the IPM model are predicted to occur at high energy and momentum. As 
a consequence, when trying to access these correlation effects with the aid of 
(e,e'p) measurements, one has to probe high missing energies, where the exclu- 
sive nature of the reaction can no longer be guaranteed. In this missing-energy 
range, the usefulness of the factorized cross section (60) becomes questionable. 
Indeed, contributions from genuine two-body photoabsorption, like e.g. those 
mediated by pion exchange, to the (e,e'p) cross section cannot be related to 
the single-particle spectral function, even when adopting a plane wave descrip- 
tion for the escaping particles. On the other hand, one-body photoabsorption 
with subsequent emission of two nucleons is a signature of ground-state corre- 
lations and has also been predicted to be the major mechanism [8] feeding the 
high-momentum components in P{k,E). Or put in other words, just as for 
the (e,e'NN) reaction, meson-exchange and isobaric currents are considered as 
unwanted noise when embarking on studies that aim at probing ground-state 
correlations with the semi-exclusive (e,e'p) and (7,p) reaction. 

In order to study the q-dependence of the semi-exclusive proton knockout 
channel we have gathered in Fig. A. 14 a number of ^^C spectra that have been 
obtained for an energy transfer of about 200 MeV. The data include both real 
and virtual photon results and are obtained at several labs. A striking fea- 
ture of the data is the ratio of the continuum strength (E^ > 30 MeV) to the 
strength residing in the discrete part of the spectrum. For the real photon case 
the continuum strength dominates the spectrum. For the highest momentum 
transfer considered here (q=585 MeV) the strength for one-nucleon removal 
from the Ip and Is shell is clearly showing up. The real photon data do not 
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show any obvious indication for one-proton emission from the Is-shell and the 
strength rises clearly beyond the two-particle emission threshold. An interest- 
ing case is the virtual photon spectrum that was obtained at low momentum 
transfer (q=270 MeV). Here the strength is a rather flat distribution that 
extends into the region of s-shell removal without reflecting a clear bump. 
Also shown in Fig. A. 14 is the calculated contribution from proton-proton 
and proton-neutron emission to the semi-exclusive spectrum. These calcula- 
tions have been performed according to the formalism outlined in Sect. 2. 2 and 
include the meson-exchange and isobaric currents in addition to the strength 
generated by the Jastrow correlations. The latter are accounted for in the 
FHNC parametrization. Apart from the highest momentum transfer consid- 
ered here, the calculations give a reasonable account of the measured strength. 
The major contribution to the semi-exclusive spectrum is ascribed to proton- 
neutron knockout. This channel is predominantly fed by the meson-exchange 
and isobaric currents, which imposes some serious limitations on the suitabil- 
ity of the semi-exclusive reaction to probe ground-state correlations. Referring 
to Fig. A. 14, it should be noted that the data taken at q=270 MeV are probing 
the dip region, whereas at q=585 MeV one is facing quasi-elastic conditions. 
One could thus conclude that semi-exclusive {'j^*\p) processes are reasonably 
well understood at the real photon point and in the dip region. The situation 
is however quite different in the quasi-elastic regime for which the origin of the 
continuum strength cannot be explained in terms of two-nucleon knockout. 

For all four situations considered in Fig. A. 14 only situation (c) exhibits some 
sensitivity to the Jastrow correlations. For the other three cases photoabsorp- 
tion on the meson-exchange and A currents dominates the calculated two- 
nucleon knockout strength. Under the kinematical conditions of (c) data are 
available for proton escaping angles varying from 27 ° up to 162 °. In Ref. [50], 
we have compared our (e,e'p) model predictions for the full range of proton 
angles. A selective sensitivity to dynamical SRC was observed. The qualita- 
tive behaviour of this sensitivity could be explained within the context of the 
two-nucleon correlation model [8] that predicts an increased sensitivity to SRC 
effects when the following relation between the missing energy and momentum 
is (approximately) obeyed 



where (E^'^ \ is the average excitation energy of the A — 2 system if the 
nucleons are escaping from the orbits characterized by h and h'. Apart from 
some factors that are related to recoil effects, the above relation reflects a 
picture in which correlated nucleons occur in pairs with respective momenta 
Pm and —pm- The kinematical conditions of the four situations considered in 
Fig. A. 14 are shown in an [Em,Pm) graph in Fig. A. 15. The shaded region 
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is the area for which the above prescription (61) would predict an increased 
hkehhood to detect ground-state correlation effects. The width of the shaded 
region was obtained by assuming that < (E^f^ ) < 50 MeV ^ thus covering 
knockout from the (Ip)^, (lp)(ls) and (Is)^ orbitals in ^^C. It was further 
assumed that correlation effects only occur for pm > 300 MeV/c. Despite 
the fact that some likelihood to observe SRC effects is predicted for the real- 
photon cases (a) and (b), the transverse character of the reaction makes the 
two-body currents dominating. For the situation (c) a proper balance between 
the longitudinal/transverse degrees of freedom and the condition (61) seems 
to be reached. 

A possible way of learning more about the short-range correlations is the sep- 
aration of the different (e,e'p) structure functions. In particular, the longitudi- 
nal response function opens good perspectives in that respect. As explained in 
Sect. 2.3, in lowest relativistic order the pion exchange effects are not affecting 
the longitudinal channel. For finite nuclei, only one experiment that separates 
the different structure functions in the high missing-energy region has been 
reported. In Ref. [54] the ^^C(e,e'p) structure functions have been separated 
in parallel kinematics {pp || q). Parallel kinematics is quite favourable for these 
purposes as only two structure functions {Wt and Wl) are contributing. The 
data of Ref. [54] are shown in Fig. A. 16 and are characterized by an excess 
strength in the transverse response function in the missing energy region above 
the s-shell peak. The transverse {Rt) and longitudinal {Rl) structure func- 
tions of Fig. A. 16 are obtained by dividing the longitudinal (transverse) part 
of the (e,e'p) cross section by (TmVc (o'mVt)- We have investigated in how far 
the observed excess strength in the transverse channel can be attributed to 
photoabsorption on correlated nucleon pairs with emission of two nucleons. In 
these calculations we have included the absorption mechanisms related to the 
SRC, IC and MEC. As becomes clear from Fig. A. 16, the present model can- 
not account for the excess transverse strength. Remark that in the presented 
model calculations all two-nucleon emission strength in the longitudinal chan- 
nel is attributed to short-range effects. 

As we find that two-nucleon knockout cannot explain the excess transverse 
strength at high missing energies, we have investigated in how far single- 
nucleon knockout from the deep-lying shells could be at the origin of the 
different missing-energy behaviour in the longitudinal and transverse channel. 
For that purpose we have calculated the cross section for knockout from the 
lsi/2 orbit at various values of Tp, covering the whole missing-energy range of 
the data. These calculations were done in a direct knockout model including 
both one- and two-body currents in the initial photoabsorption process [55]. 
For the whole missing-energy range covered here {Em < 70 MeV) we use the 
bound state lsi/2 wave function as obtained from a HF calculation. Formally, 
this means that in the calculations the spectral function 
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n 

x6 {e - {Ef - E^)) , (62) 

is used in a factorized form : 

Pi,{k,E) =1 (Pi,{k) p 5,,(i?)(2j + 1) . (63) 



The function Sij describes the spreading of the hole strength characterized by 
the quantum numbers Ij as a function of the missing energy in the A-1 system 

SiAE) = T.\{^n-\E§-'')\ci,\ME^))"6(E-(Et'-E^)) , (64) 



and (t)ij{k) is the mean-field single-particle wave function in momentum space. 
Microscopic calculations of the single-particle spectral function [9], however, 
show that deviations from mean-field wavefunctions predominantly occur at 
higher missing energies {Em > 50 MeV). As such, the factorized prescription 
(64) could be considered as a reasonable approximation as long as relatively 
low missing energies are probed. In order to account for the spreading of the 
Is hole strength in ^^B we rely on a prescription which is due to Jeukenne 
and Mahaux [56]. In this parametrization the nucleon hole spectral function 
is determined by a Lorentzian with an energy-dependent width T[E) : 

g i^\ ^ i ZijT{E) 



where Zij (eij) is the quasi-particle strength (energy) of the hole state under 
consideration. The width T{E) is then related to the imaginary part W of the 
optical potential (F = 2W) for which the following parametrization is adopted 



9.iE-S, 



N 



\4 



With the above assumptions a fair description of the missing-energy depen- 
dence is reached for the longitudinal (e,e'p) results shown in Fig. A. 16. It is 
worth mentioning that the fact that the Is strength has a high missing energy 
tail can be attributed to correlation effects. In order to reach the agreement 
of Fig. A. 16 a spectroscopic factor of 0.5 for the lsl/2 state has been adopted. 
Despite the reasonable description of the longitudinal structure function, the 
transverse strength remains underestimated over the whole missing energy 
range. Final-state interactions are a possible explanation for this observation. 
Eventhough differences in the way FSI affects the different response functions 
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have been noted, the FSI effects are unhkely to act in such a selective manner 
that they produce a long tail of strength in the transverse response function 
leaving the longitudinal one almost unaffected. 

In Fig. A. 17 we show a similar theoretical analysis for the measurements of 
Ref. [53] that were performed at higher energy and momentum transfer. The 
^^C(e,e'p) data of Fig. A. 17 refer to quasi-elastic conditions and were obtained 
in parallel kinematics (pp||^. As the data cover a larger missing-energy range 
than those of Ref. [54] the approximation (63) is at stake. The data show 
a clear bump in the missing-energy region that is dominated by proton re- 
moval out of the Is-shell. Above the s-shell region, the data exhibit a long 
tail. Part of this strength can be attributed to two-nucleon knockout. For 
the largest momentum transfer considered here, however, the calculated two- 
nucleon knockout strength is just a small fraction of the measured strength. At 
q=585 MeV/c a reasonable account of the data is achieved when adding the 
calculated 2N knockout strength to the contribution related to knockout from 
the Is shell. In the region just above the Is bump where the approximation 
(63) is expected to be reasonable, the calculations underestimate the data in 
both cases. It remains to be investigated whether model calculations starting 
from realistic spectral functions can explain the origin of this strength. 



4 Conclusion 



In this paper we have provided a framework in which the effect of ground- 
state correlations on electroinduced one- and two-nucleon knockout cross sec- 
tions could be quantified. We have gone beyond the standard shell-model 
by considering short-range corrections to the nuclear wave functions. In the 
first-generation calculations presented here, this has been achieved by a state- 
independent Jastrow ansatz. The effect of the short-range effects on the cross 
sections is estimated on the basis of a cluster expansion. It is shown that when 
considering correlated wave functions, standard one-body photoabsorption in 
the Impulse Approximation generates a whole chain of "effective" absorption 
mechanisms that are 2,3,..., A-body in nature. These multi-body effects can be 
related to the ground-state correlations and compete with regular two-nucleon 
photoabsorption mechanisms, as e.g. generated through meson-exchange and 
A33 creation. 

In the model outlined, a partial-wave expansion technique is adopted for 
the description of the final-state interaction in A(e,e'NN) processes. A di- 
rect knockout reaction model is adopted. Special care is taken to treat all 
spin degrees of freedom and anti-symmetrization effects exactly. The model is 
particularly suited to calculate exclusive A(e,e'NN) angular cross sections for 
creation of the residual nucleus in a specific state. 
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A novel technique to calculate the 2N knockout contribution to the semi- 
exclusive (e,e'N) channel is outlined. With the aid of Racah algebra, the in- 
tegration over the solid angle of the undetected escaping nucleon could be 
done analytically. This procedure cuts severely on the amount of numerical 
integrations which have to be performed. The procedure of calculating the 
semi-exclusive (e,e'p) strength adopted here, has the advantage of being unfac- 
torized in nature. This allows calculating the contribution from ground-state 
correlations without excluding contributions from other sources, as e.g. these 
arising from photoabsorption on two-body currents. Such contributions are 
completely excluded when expressing the semi-exclusive (e,e'p) cross section 
in terms of the one-body spectral function. 

As far as the exclusive 2N knockout channel is concerned, the numerical re- 
sults presented here include both (e,e'pp) and (e,e'pn) angular cross sections 
for the target nuclei ^^C and ^^O. The shapes and magnitudes of the angu- 
lar cross sections reflect the complexity of the photoabsorption mechanism 
and the shell-model structure of the target and residual nucleus. The (e,e'pn) 
channel is characterized by considerably larger cross sections than the (e,e'pp) 
channel. Electroinduced proton-neutron emission is dominated by the mesonic 
and A-isobar degrees of freedom, with the central ground-state correlations 
playing only a marginal role. Therefore, the proton-neutron knockout strength 
is predominantly transverse in nature even when longitudinal kinematics is 
adopted. Two-proton knockout, on the other hand, exhibits a selective sen- 
sitivity to ground-state correlation effects. A major source of electroinduced 
two-proton knockout is, however, the AA^ -^ pp mechanism. It is therefore 
essential to have a realistic description of these isobar degrees of freedom. 
This requires the input of knowledge obtained within the context of pion and 
real-photon absorption on nuclei. 

The conclusions drawn for the (e,e'NN) reaction have their implications for 
the semi-exclusive (e,e'p) channel as correlation effects will manifest them- 
selves predominantly as two-nucleon knockout. The calculations predict this 
channel to be dominated by proton- neutron knockout, which in its turn is 
strongly fed through the meson-exchange and isobaric currents. Apart from 
the two-body currents considered here, also multi-scattering effects can con- 
taminate the link between the semi-exclusive (e,e'p) data and the spectral 
function. This effect has been the object of several investigations [57,58] with 
sometimes different conclusions. All this puts heavy constraints on the ap- 
plicability of the semi-exclusive (e,e'p) reaction to gain empirical information 
about the single-particle spectral function. In any case, the strength generated 
by meson-exchange and A degrees of freedom has to be carefully estimated 
and substracted from the (e,e'p) data before relating the measured strength 
to ground-state correlation effects. Better conditions to detect the correlation 
effects in (e,e'p) are predicted to occur at higher momentum transfer and 
specific bands in the {Em,Pm) configuration space. 
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A Two-body matrix elements with short-range currents 



The purpose of this Appendix is to give the expressions for the two-body 
matrix elements (Eqs. (21) and (20)) with the effective two-body currents 
that account for the short-range effects in the initial wave function (Eqs. (53) 
and (54)). As outlined in Sect. 2.4 these effective operators are composed 
of the one-body operators of the impulse approximation and a Jastrow like 
correlation function that corrects for short-range correlations in the initial- 
state wave function. 

The longitudinal strength attributed to SRC is contained in the two-body 
charge density of Eq. (54) and is determined by the following effective Coulomb 
operator (17) : 



MTm (p|L(1,2)) =- (jj(gri)r7M(fii)^(ri2)e 
+ 3j{qr2)yjM{^2)g{ri2)e 



1 + r,,i 
2 

1 + r. 2 



(A.l) 



With the aid of the expansion (55) this leads to the following two-body matrix 
elements : 



<a6;Ji||Mr'(p|t^(l,2))||crf;J2>= 

- J2 V^e^^ <10L0\J0> dn dr^g'iir-i, r^) 

IL ^ J J 

X [5ac,iT X{ja, jb, Ji;jc, jd, -h] L, /, J)D{a, c, L, ri)D{b, d, I, r2)jj{qri) 
+Sbd,7T X{ja,jb, Jujcjd, -h] I, L, J)D{b, d, L, r2)D{a, c, /, ri)jj(gr2)] (A.2) 



with 



D{a,b,l,r) =< laja\\Yl{n)\\lbjb >r 



(A.3) 



The radial transition density < a||0||6 >r is defined such that it is related to 
the full matrix element through < a||0||6 >= J dr < a\\0\\b >j.. In the above 
expressions X denotes the 9j symbol in the conventions of ref. 



X{ja, jb, Jl] jc, jd, J2; jejf, J) = < 



Introducing the operator 



ja jb Jl 

jc id J2 

Je Jf J 



> . 
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0}^(g)= Y. fdf<J + KMilM2\JM> 

Mi,M2 

xFj+«Mi(fi)jj+«(gr)JM2(r) , (A.4) 

the electric and magnetic transition operators (Eq. 16) can be rewritten as : 

TTM%<l) = OrM' (A.5) 



i(-l)^-'+ 



TfMiQ) = E ' 7 ' ^J + K-iO^JM ■ (A.6) 

Accordingly, the matrix elements of the operator O}^ suffice to determine 
both the electric and the magnetic strength. The transverse strength at- 
tributed to SRC is determined by the effective two-body current of Eq. (53). 
The latter has two components related to the convection and magnetization 
one-body current of the impulse approximation. For the component related to 
the magnetization current the effective O'jj^ operator reads : 



0} (^'ir'"!!- 2)) = - / rff ^ < J + /€ Ml 1 M2 I JM > 

•^ M1M2 

^ {^^'^(^"^1)^(^12) [(y X ai) ^Yj+^MA^)3J+n{qryj 

+ ^^i^- ^2)9{ru) ((V X a,)^^ lj+.M.(fi)jj+«(gr))} (A.7) 

After some lenghty but straightforward manipulations, the transition matrix 
element for this operator can be reduced to : 

< ab; Ji \\0^j (4f "(1> 2)) \\cd; -h >= v^ '""^ 



Mm 

J\IjJ2J J A 



^ J2 J2 ^ — - <lOLO\J + K + riO> 

J 1 J + K + ri\ ( J + K J + K + Tj 1 



r)=±l lL.h ^ 



x^J + «: + Vi^^ / J4 i I 1 1 J 

X J dn J dr2gi{n,r2) {[fip5ac,n + f^Jac^K-lf^'^^-^^ 

xX{ja,jb, Jujcjd, J2] Ja, I, J)B{a, c,J + K + r],L, J4, qri)D{b, d, I, ra) 

+ [fl'pSbd,n + fJ'nSbd,u] ^ija, jb, Jl'ijc, jd, J2'i I, Ji, J) 

xB{b,d, J + K + T], L, J4,qr2)]D{a,c,l,ri)} , (A. 

where we have introduced the radial transition density 

B{a, b, I, Ji, J2,pr) = < a\\ji{pr) [Yj^Q) ® a]j^ \\b >r 
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= jaJ2JbX{la, 1/2, ja] k, 1/2, jb] ^1, 1, ^2)/a^lJ — 

xi-lY' <laOJi 0\h > r^Mr)Mr)ji{pr) 

The matrix element for the SRC contribution related to the convection current 
can be derived in an analoguous manner and reads : 



K ( l{2],conv 



<ah-.h\\0^j[.Jfj:^-\l,2))\\cd-.h>= 

I — 6 X — ^ Jt J'2J JaJ-/ , , 

V^— — y \. </OLOJ + /tO> 

^Mn^^, I 

X {6ac,7rX{ja,jb, J^jcjd, J2] J4, 1, J)E{a, c, J + K, L, J4, qri)D{h, d, I, r2) 

+<5m,.(-1)'+^*+^X(j„ jfe, Ji; je, Jd, J2; /, J4, J) 

xE{h,d,J + K,L,Ji,qr2)D{a,c,l,ri)} , (A. 9) 

where we have introduced the radial transition density E{a, b, I, Ji, J2,pr) 

E{a,b,l,Ji,J2,pr)^<a\\ji{pr)\Yj,{n)(g){V-V')]^ ||6 >, . (A.IO) 



The notation V' refers to a gradient operator acting to the left. The terms 
involving the derivates of the central correlation functions have been neglected 
in the above expression. 
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Fig. A.l. Adopted conventions for the angular variables of an A(e,e'N(jN6) reaction. 
The electron scattering plane is the xz plane. 




Fig. A. 2. Two-nucleon knockout as interpreted in a shell-model picture. After ejec- 
tion of two nucleons with energy £„ and eb the residual nucleus is left in a two-hole 
state (hh')"^ 
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Fig. A. 3. The four diagrams corresponding with A33 creation and vr-exchange. 
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Fig. A. 4. Different contributions to photoinduced two-nucleon knockout. Meson 
exchange processes are denoted with the dashed line. The photon is indicated with 
a wavy hne and the ground-state correlations with the heavy dots. 
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Fig. A. 5. Differential cross sections for the ^^C(e,e'pp) reaction at e=700 MeV, 
cj=250 MeV and 9e=30° (q=383 MeV/c) for excitation of an | (lp3/2)"^;0+ > and 
(lP3/2)~^;2+ > state respectively. 
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Fig. A. 6. As in Fig. A. 5 but now for the C(e,e'pn) reaction 
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Fig. A.7. Differential 

cross sections for the "'^^0(e,e'pp) reaction at e=1.2 GeV, ix'=225 MeV and ^e=12° 
(q=319 MeV/c) for excitation of an | (1^3/2) ""'^(1^1/2)""'^; Jr = "^^ > two-hole state. 
(a) full cross section (b) longitudinal contribution. 
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Fig. A. 8. As in Fig. A. 7 but now for proton-neutron emission. 
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Fig. A. 9. Differential cross sections for the ^^0(e,e'pp) 
reaction at e=1.2 GeV, tj=375 MeV and 6le=12° (q=429 MeV/c) for the excita- 
tion of an I (1^3/2)"^ (lPi/2)"'^; <^R = 2"*" > two-hole state, (a) full cross section (b) 
longitudinal contribution. 
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Fig. A. 10. As in Fig. A. 9 but now for proton-neutron emission. 
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Fig. A. 11. (a) Coplanar and symmetrical kinematics, (b) the c.o.m. momentum P 
as a function of the opening angle 6 for two-proton knockout from ^^O at an energy 
transfer of u;=210 MeV. 



46 



0.175 - 



0.15 



0.125 



■a 



a 

T3 

a 

a 
a 

D 

T3 



0.1 



0.075 



0.05 



0.025 - 



^'^0(e,e'pp)^^C 



(lp3/2)-' 



(lp3/2)'(lpl/2)' 



(ipi/2r 




ep(deg) 



Qp (deg) 



ep(deg) 



Fig. A. 12. The ^^0(e,e'pp) differential cross section as a function of the proton angle 
in coplanar and symmetrical kinematics for e=516 MeV, uj=210 MeV and ^£=31". 
The different combinations for emission out of the p-shell orbits are considered. 
Different central correlation functions have been used : OMY (dotted line), FHNC 
(solid line) and MC (dashed line). The dot-dashed line is the result of a calculation 
in which only the isobaric currents are included. 
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Fig. A. 13. The four terms contributing to the ^^0(e,e'pp)((lp3/2)^"'^(lpl/2)^^) an- 
gular cross sections plotted in the central panel of Fig. A. 12. The dot-dashed line 
shows the result when including only the delta-currents. The other curves are ob- 
tained after adding also the central short-range effects. Different central correlation 
functions have been used : OMY (dotted line), FHNC (solid line) and MC (dashed 
line). 
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Fig. A. 14. Missing-energy dependence of the semi-exclusive ^^C(e,e'p) and ^^C(7,p) 
spectrum for u ~ 200 MeV and various kinematical conditions. The dotted (dashed) 
Une shows the calculated contribution from proton-proton (proton-neutron) emis- 
sion. The solid line is the incoherent sum of all the calculated strength contributions. 
The dot-dashed line is the calculated contribution from two-nucleon knockout in- 
cluding only the transverse two-body currents (for most cases this line falls on top of 
the solid one). The data are from Refs. [14] (circles), [53] (squares), [51] (triangles) 
and [52] (stars) . 
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Fig. A. 15. The missing nionientum as a function of the missing energy for the 
kinematical conditions of Fig. A. 14. 
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Fig. A. 16. The longitudinal and transverse structure functions versus missing energy 
for the ^2c(e,e'p) reaction at uj=l22.5 MeV and q=397 MeV/c. The dashed line is 
the calculated contribution from proton-proton and proton-neutron emission. The 
solid line is the predicted contribution for one-proton knockout from the Is shell 
applying the same spectroscopic factor in the longitudinal and transverse structure 
function. The data are from Ref. [54]. 



51 




C(e,e'p) (0=210 MeV q=585 MeV 

(a) 



20 40 60 80 100 120 140 

missing energy (MeV) 




50 100 150 200 

missing energy (MeV) 



Fig. A. 17. Missing-energy spectra for ""^^0(6, e'p) at quasi-elastic kinematics, (a) 
e=505.4 MeV, uj=2W MeV and q=585 MeV/c (b) e=698 MeV, a;=355 MeV and 
q=775 MeV/c. The dotted (dashed) hne is the calculated strength from (e,e'pp) 
((e,e'pn)). The dot-dashed line the calculated strength from exclusive one-proton 
ejection out of the Is shell. The solid is the summed strength from all calculated 
contributions. The data are from Ref. [53]. 
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